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1 Introduction and summary of results 



In its original and most studied form p], the gauge-string correspondence relates a field the- 
ory where all fields transform in the adjoint representation of the gauge group to a theory of 
closed strings. An enormous amount of work has been devoted to generalize the correspon- 
dence in different directions. Since the physics of real world strong interactions is dominated 
by mesons and baryons, a generalization of obvious importance is to include matter trans- 
forming in the fundamental representation of the gauge group. Early works on the subject 
were 0121111. 

In the seminal papers [3 IH] , it was argued that including a small number of flavors {Nf 
fixed, Nc — > oo) corresponds in the dual theory to adding an open string sector due to the 
presence of a few flavor branes. In the field theory, this is a quenched approximation in the 
sense that the diagrams with fundamentals running in the loops are suppressed by N~^. In 
the gravity side, one can neglect the backreaction of the branes on the geometry. These ideas 
and methods sparked a lot of activity which led to a neat understanding of many different 
phenomena in such a limit, in a variety of setups. 

Nevertheless, field theories in which Nf is of the same order as Nc are amazingly rich as is 
clear, for instance, from Seiberg's analysis of AT = 1 SQCD. In addition, the fact that in the 
real world ^ is of order one is reflected in experimental facts as, among others, multihadron 
production or the large mass of the rj' meson. This motivates the study of string duals in the 
Veneziano limit Nc ^ oo with g'^Nc and ^ fixed. In order to maintain the curvature small, 
strong coupling (large g'^Nc) is required. However, this case is not as well understood as the 
quenched one. Some works along this direction studying four-dimensional field theories are 
121 ini El ISl Uni UH 1121 (three-dlmenslonal field theories were addressed in [121 CH US] ) . In one 
way or another, all these models have a region of large curvature associated to the presence 
of the fundamental quarks. 

In this note, a setup consisting of wrapped branes, dual to A/" = 2 SQCD is presented. By 
appropriately smearing the flavor branes such that the shell of branes forms a domain wall 
in the geometry, they do not generate any pathological region in the geometry (as long as 
Nf < 2Nc so there is no Landau pole in the dual field theory). This smearing corresponds to 
a smearing of the complex masses of the hypermultiplets, more details are given in section 
13.21 Anyhow, the solution has a naked curvature singularity in the IR and a divergent dilaton 
in the UV. These features are already present in the original setup and are not related to 
the presence of flavor branes. Some comments on how these pathologies can be dealt with 
can be found in sections |21 and jH respectively. 

The framework is the gravity dual of A/" = 2 SYM introduced in J^l Ej ■ It consists of D5- 
branes wrapping an inside a CY2 in such a way that eight supercharges are preserved. At 
small energies compared to the inverse radius of the sphere, the theory becomes effectively 
four- dimensional. It will be shown how one can add a different kind of D5-branes in a 
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way that no further supersymmetry is broken. They introduce an open string sector that 
corresponds to having in the field theory (massive) hyper mult iplets transforming in the 
fundamental representation of the gauge group. As explained above, when the number of 
these branes is comparable to the number of colors, their backreaction on the geometry has 
to be taken into account. This is done by applying the procedure advocated in [12]: the D5- 
branes are sources for the Ramond-Ramond F(3) field strength, and therefore they modify 
its Bianchi identity. Once this is taken into account, the set of first order equations can 
be found by imposing the vanishing of the supersymmetry variation of the fermion fields. 
The resulting system of equations is quite complicated and I have not managed to find an 
analytical solution. However, as it will be shown, it can indeed be explicitly solved in a 
certain region of the space. This is enough to prove that the modification of the couplings 
due to the fiavor branes is the one expected from field theory. Finally, the meson spectrum 
is studied. The excitations can be rearranged in a tower of massive Af = 2 vector multiplets. 
Due to the UV pathological behavior of the solution, a regularization procedure is needed in 
order to estimate the masses of the physical excitations. Once this procedure is implemented, 
the dependence of the meson masses on the different physical parameters is briefly discussed. 

Apart from the interest of the particular model that will be addressed, one of the aims of 
this paper is to make clearer the methods of ^2] by their application to a simpler setup. In 
particular, as stated above, in this Af = 2 construction the flavor branes do not spark the 
divergence or vanishing of any component of the metric. Hopefully, this note can contribute 
in the development of techniques for the study of holographic duals with unquenched matter. 

The contents of the paper are organized as follows: in section |21 the unflavored solution is 
reviewed and rewritten with a notation that will be useful in the following. In section IH.H 
the embedding of the flavor branes is discussed and the physical meaning of the different 
spacetime directions is explained. In section 13.21 the system of equations including the 
backreaction is obtained and in section EiHl it is shown how the couplings are modifled by the 
unquenched flavor. Section |3] deals with the spectrum of excitations. A few WKB formulae 
used in section are summarized in appendix 1X1 

2 A review of the dual of A/" = 2 SYM 

In references fHl EI 5 a HB gravity dual of A/" = 2 SYM in the Coulomb branch was found. 
The ten-dimensional solution was obtained as an uplift from seven- dimensional 5*0(4) gauged 
supergravity. It corresponds to 5-branes^ wrapping a two-sphere with the appropriate twist- 
ing to preserve eight supercharges, i.e. A/" = 2 in the effective four- dimensional low energy 
theory. Geometrically, it corresponds to wrapping the branes along a compact SLag two- 

^ Originally, the solution was built with NS5-branes. In the following, the S-dual solution corresponding 
to D5-branes will be considered. 
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cycle inside a Calabi-Yau two-fold. This leaves two fiat transverse dimensions which are 
identified with the moduli space corresponding to giving vevs to the complex scalar inside 
the M = 2 vector multiplet. 

The solution presented below has a curvature singularity in the IR. The singularity is good 
according to the criterion of In fact, it was shown in jTHI, by considering wrapped NS5- 
branes that the singularity is an artifact of the supergravity approximation and is resolved 
by the worldsheet CFT. Therefore, the singularity will not be a matter of concern in the 
following. 

This section does not contain original material but reviews the results of [THt IT7| ED] in 
order to fix notation for the following and make this note reasonably self-contained. An 
excellent review of this model and similar constructions is 1211. 



2.1 The solution 

The metric (Einstein frame) of the solution found in EI reads: 

1 



NcQsOi e 



2 



-dxls + z{de' + sin' Od^') + e'^'dz' + de'+ 



+ 



cos 



1 + COS 9d(p) + — sin 



(2.1) 



there is a magnetic RR twoform: 

C(2) = Ncgsa'(j)2d 
and the dilaton is given by: 



sin^ e 



— {d(f)i + cos 6d(p) 



(2.2) 



sin^^ 



1 + ce" 
27 



-2z 



(2.3) 



and we have defined the quantities: 



-2x 



1 - 



1 + ce 



-2z 



2z 



n = e'' cos' 9 + e"^ sin^ 9 



(2.4) 



The angles take values in the following intervals: 9 G [0,7r], 9 G [0, |], 0i,02 £ [0,27r). 
The variable z ranges from a value zq for which e~'^^^°^ = up to infinity. 

There are two integration constants on which the solution depends, $o and c. Following 
the discussion of [121 Ej, the solution is dual to A/" = 2 SYM at points of the Coulomb branch 
in which the vevs for the entries of the complex scalar matrix are distributed in a ring (they 
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have fixed modulus whereas the phase is smeared). The constant c > — 1 determines the 
size of the ring. Solutions with c < — 1 are unphysical and, in fact, the singularity becomes 
of the bad type. Notice that although taking a single vev different from zero spontaneously 
breaks the U{1)r, this smearing procedure restores it, at least if one ignores A''"^ effects, 
where one could start seeing that the distribution is not continuous (see figure [T] in section 
13. 2|) . Points of the Coulomb branch where the U{1)r is not restored were discussed in [T7j. 
The quantization condition is: 

^ / ^(3) = N.T, (2.5) 
^^(10) Js^ 

where is the transverse three-sphere parameterized by ^, 0i, 02 at 2; — oo, yielding: 

F^3) = Ncgsa' d^i d^2 2 sin 9 cos 9d9 = in^N^g^a' (2.6) 
53 Jo Jo Jo 



Taking into account: 



1 11 

= {2n)%a'^ ' 24;;; = (2^oW ^^'^^ 



one can immediately check that ()2.5p is satisfied. 



2.2 Rewriting the solution and Killing spinors 



In reference the solution was rewritten in different variables which allow a better un- 
derstanding of the physics. Define: 



p = sm I 



a 



zcosl 



(2. 



so the metric ()2.H) reads: 

dsjf) = QsNca'e 



-dxf ^ + z{d9^ + sin^ 9d(p^) + 



—2^ 

+e-^*(rfp^ + p^dcpl) + [da^ + ^^(#1 + cos 9d<ff^ 



(2.9) 



Written in this way, it is clear that the Calabi-Yau twofold directions are 9, ip, a, (pi (of course, 
in this solution with fluxes there is not a Calabi-Yau any more, but it can be thought of 
as a deformation of the Calabi-Yau that was present before backreaction) . The coordinates 
p, 02 span the transverse two-dimensional plane, so they should be identified with the moduli 
space, and therefore rotations in 02 are related to the U{l)Ft symmetry of the field theory. 
These statements will be made more precise below. 
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Actually, one could start by writing the metric ()2.9|) as an ansatz with z{p,a), $(p, cr). 
The ansatz for the RR 3- form would be: 

F(3) = NcQsO.' —g'd(t>2 A dp A {dcpi + cos 9dip) — gd(f)2 Ada A {dcpi + cos 6d(p) + 

+gsmed(f)2 Ade Ad<p (2.10) 

which ensures (iF(3) = and we have defined: 

' = dp, ■ = d^ (2.11) 
Introducing this ansatz in the IIB transformations for the fermions: 

= a^e + ^^r'^^e + ^et (r™-95;:^r^^'^3)^^,^^^^^^^ (2.12) 

a set of first order equations can be obtained. Consider the orthonormal frame: 

= e^dxi {i = 0, ...,3) , = \fgsN^' e'^ ^dQ , = ^/ QgNcO'e^ \/z sinOdcp , 
e = \/gsNca'e * dp , e' = ^Jg^N^a'c * p(i02 , 

= ^/g^K^'^-^do , = >/gJ^^—Lcr{d4)i + cosdd(p) (2.13) 



The Killing spinors are given by: 

e = ete"^'^'^'*"+^'^'^"77 (2.14) 
where is a constant spinor satisfying: 

(yJ^ib&lV = -V , ^45891] = v (2.15) 

The first order equations read: 

g = -pz', (2.16) 
a 



e^* = — , (2.17) 

zz 



2$ 

g' = -26-2* pa<i> , (2.18) 
g = -z-^e~^^ag + 2z-^pae-^^^' . (2.19) 

It is possible to check that these equations ensure the equation of motion for the 3-form 
(i(e* * F(3)) = 0. Notice, however, that these equations are not independent since fj2.19|) is 
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automatically implied by ()2.16p and ()2.17p . The system ()2.16|) - ()2.19|) can be rephrased as a 
single second order non-linear partial differential equation (PDE) for the function z{p, a): 

pz{z - oz) = a{pz^ + z' + pz") (2.20) 

and once z is known, g and $ can be obtained from ()2.16|) . ()2.17|) . The implicit relation 
deduced from ()2.8j) : 

a2-ze-2"(e2"-p2) = (2.21) 

(where e~^^ was defined in fl2.4|) ) solves ()2.20|) and yields the known solution'^ ()2.H) - ()2.4|) . 
The value of g in terms of the old variables can be read by comparing ()2.2|) to ()2.10|) : 

sin^ 6 

qZx pQgZ fj _|_ gjj^ y 

It is also useful to define quantities zq, po which depend on c as: 

e-2x(.o) _ , Po = (2.23) 

The minimum value that z can attain is zq- When c = — 1, then zq = 0. 

Using ()2.9jl and ()2.15|1 . it is immediate to see that two kinds of D5-brane probes can be 
added to the setup preserving the full supersymmetry of the background. They correspond 
to the K-symmetry projectors O"iroi2345 and criroi2389- The first kind corresponds to D5 
extended along xq, . . . , xs, 9, if placed at a = (notice that the dilaton does not diverge at 
a = despite the appearance of ()2.17j) since i ~ a near cr = as can be deduced from 
fl2.20|) ). These branes, which wrap a compact S"^, are the color branes. Computing the 
Born-Infeld action on their worldvolume leads, upon reduction to four dimensions, to the 
J\f = 2 SYM action. In |2i]|, the identification of the Yang- Mills coupling, the 9ym and the 
radius energy relation was made precise. In order to follow their argument, notice that at 
a = 0: 

{a = 0, p<po)^r^l' , (a = 0, p>po)^r^^°f^ , (2.24) 

where we have used (j2I2H), ()2.1(ij) in order to derive these expressions, cr = 0, p < po is the 
singular locus of the geometry. As acknowledged in jT^, it is remarkable that one obtains 
regular results for physical quantities despite the singularity. 

Now, consider the abelian worldvolume DBI-I-WZ action of a D5-brane probe placed at 
(7 = 0. Then, expand it to leading order and promote it to a non-abelian one. Normalizing 
the SU{Nc) generators as tr{T^T^) = ^f^^'^, one obtains 



Sym = / d'x{\F^,Ff + ^D^^^D^^a} + / d'xF^.Pf (2.25) 



■^Simple solutions of 12.20|l are z ~ ci + C2 logp which is not physical since from (|2.17() one gets $ = oo 
and z — c\\j2(j^ + C2 which leads to constant dilaton and g — Fi^y^ — and just Ricci flat Ri_5 x EH^ where 
EHa denotes the Eguchi-Hanson space. 
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where we have defined which is the complex scalar of the M = 2 vector multiplet as: 



-pe'^^ (2.26) 

After integrating the S"^ parameterized by 9, (p, the value of the couplings can be read from 
the solution ()2.24j) . For p < po they are constant while for p > po they read: 

'(^U=o) = T4logp, (2.27) 



Oym = 2((7U=o)Ar,02 = -2Ar,02 (2.28) 

Now, noticing that the complex scalar \E' has protected mass dimension one, one is lead, from 
()2.26|) to the following radius-energy relation: 

p = ^ (2.29) 

where p is the mass scale at which the theory is defined and A the dynamically generated 
mass scale. Then: 

(i[9YM) = 01 / /AX = -w] = -^9ym 2.30) 

c'log(/i/A) alogp Svr^ 

On the other hand, one can compute the chiral anomaly, which is associated to U{1)r shifts: 

02 ^ 02 + 2e Oym ^ Oym - 4:N,e , e G [0, 27r) (2.31) 

The parameter e takes values in [0, 27t) since, although in the geometry one comes back to 
the same point after a rotation of 27r in 02, the complete U{1)r rotation is of a Air angle. 
This can be seen from ()2.14p . since the spinors pick a minus sign upon 02 — * 02 + 27r [16] . 
This is also consistent with the complex scalar \E' having R-charge 2, see ()2.2(ij) . 
Therefore, allowing 6ym to change by a multiple of 2n one finds the values of e: 

e=^k, k = 0,l,...,4Nc-l (2.32) 

describing the anomalous breaking U{1)r I^aNc- 

Let us finish the section by defining the holomorphic coupling r = i-l^ + which, in 

9ym -^^ 

view of (EI2Z|),(EI2H1), and defining u = Ape'^\ Uq = Apo is: 

■ Nc-. uo f I I ^ 
r = z — log — , tor m < Mo , 

71 A 

T = i — ^ log — , for \u\ > Uq (2.33) 

TT A 

which is precisely the result expected for a ring-like distribution of vevs [TB] . 



7 



3 The dual with flavor 



3.1 Flavor branes in the solution 

The second kind of supersymmetric embeddings for D5-branes is associated to the /t-symmetry 
projector (Tiroi2389- It corresponds to branes extended along xq, Xi, X2, X3, a, 0i at any point 
in the other coordinates 6, (p, p, 02- These branes provide fundamental hypermultiplets in or- 
der to build A/" = 2 SQCD. First of all, since they are extended in the a direction, their volume 
is infinite (compared to the color branes), making exactly zero the effective four- dimensional 
gauge coupling living on them, and therefore providing a global symmetry group U{Nf) if 
they are placed on top of each other. Remember that this (and not U{Nf) x U{Nf)) is the 
flavor symmetry of the Af = 2 SQCD lagrangian, due to the coupling (in A/" = 1 notation) 

Tr{Q^Q) (3.1) 

where the Q,Q represent the ( ant i) -fundamental chiral multiplets inside the M = 2 hyper- 
multiplet and ^ is the chiral multiplet in the adjoint inside the M = 2 vector multiplet. This 
is consistent with having only flavor D5-branes (anti D5-branes would break supersymme- 
try). It can be seen that the coupling (jH.lj) appears on the worldvolume theory of a probe 
color D5-brane. The easiest argument is that A/" = 2 is preserved, so (jH.lj) is needed. Also, 
from the point of view of the brane intersection, it is what one expects, since ^ parameterizes 
the two directions which are orthogonal to both color and flavor branes (those spanned by 
p, ^2). The values of p, 02 at which each flavor brane is placed give the complex mass of the 
corresponding hypermultiplet (taking into account the identifications of section IT^ . The 
position in the 6', coordinates does not play any role in the low energy field theory, since we 
are interested in energies below the inverse radius of the S'^ . However, this statement should 
be taken with certain criticism, since, as usually happens in the wrapped brane models, one 
cannot really achieve a separation of scales between the unwanted Kaluza-Klein modes and 
the relevant physical modes. A discussion on this point and some ideas of how to deal with 
this problem can be found in |22j. 

Let us now focus on the ^7(l)i? symmetry that, as described in section |21 corresponds to 
rotations of the 02 angle. Massless flavors would correspond to branes located at the origin 
of the moduli space p = 0, a point which is invariant under 02 rotations, and therefore 
preserves the U{X)r- A brane located at some p > would correspond to a massive flavor, 
and, as expected, breaks explicitly the U{X)r- 

There is an 5*0(3) isometry which acts on ^, 0i (it is not SU{2) since 0i takes values 
between and 27c whereas in the usual SU{2) left invariant one-forms, it would range up to 
47r), which does not play a role in the low energy theory Thus, the SU{2)ji symmetry 
of the field theory cannot be entirely realized in the geometry. Nevertheless, there is another 
isometry in the solution, the rotations in the 0i angle, which can be identified with the 
U{l)j C SU{2)ji. An indication that this is the case comes from comparison to a non-critical 
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string setup. It was argued in (see also [21]) that, if one considers the eight-dimensional 
background Ri^s x IR2 x ^^i^'^^'^ ^-^d places D3-branes extended along the Mi,3, at the tip of 
the cigar, and D5-branes extended along the Mi^a x — u^T)'' ^^en the low energy theory living 
on the D3 is = 2 SQCD. Heuristically, one can think of this non-critical setup as the one 
described in this paper where the S"^ parameterized by 9, (p has shrunk to string scale size. 
Naturally, the curvature would be of the order of the string scale. In |23], it was shown that, 
by considering the fermions, this U{1) forms part of an SU{2) symmetry. The wrapping of 
the flavor D5 branes does not break this symmetry, regardless the hypermultiplet is massive 
or not. However, the scalars inside the hypermultiplets are charged under the SU{2)fi, so 
the meson-like excitations lie in non-trivial representations of this group, as will be shown 
in section^ 

When two hypermultiplets have the same mass, it is possible to enter a Higgs branch by 
giving expectation values to Q, Q. From the gravity side, this should be described by using 
the non-abelian brane worldvolume action along the lines of |26j. 

Table n summarizes the setup. 



CY2 M2 





Xl,3 


cr 


01 


9 




P 


02 


D5 


















Nf D5 



















4D spacetime 




U{l)j C SU{2)r 






Energy scale 


Uil)R 




(— oo, oo) 


[0,oo) 


[0,2n) 


[0,7r] 


[0,27r) 


[0,oo) 


[0,27r) 



Table 1: A scheme of the setup: for the brane configuration, a line — means that the 
brane spans a non-compact dimension, a point • that it is point-like in that direction and a 
circle Q that it wraps a compact cycle. The physical meaning of the different dimensions 
is summarized as well as their ranges. Above, it is shown which directions spanned the 
Calabi-Yau and which the transverse plane before backreaction. 



3.2 Computing the flavor brane backreaction 

Now that the embedding for flavor branes has been identifled, the procedure of ^2] can be 
followed in order to flnd a backreacted solution and the effects of having Nf ~ A''^ in the field 
theory. The unquenched solution is a solution of supergravity coupled to the flavor branes 
which, unlike the color branes, do not disappear into fluxes. The global flavor symmetry of 
the field theory corresponds to the gauge symmetry on the brane worldvolume. As explained 
above, this symmetry becomes global, from the four- dimensional point of view because the 
flavor branes have inflnite volume since they are extended in the non-compact a direction. 
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Clearly, the solution must depend on the masses one chooses for the fundamental hyper- 
multiplets (the position in p, 02 of the branes). The case addressed below is the simplest 
one: each single mass breaks the t/(l)_R symmetry but since Nf — > oo, we may consider a 
setup in which the masses are smeared and distributed in a ring, such that the f/(l)_R is 
restored. The flavor symmetry is explicitly broken to U{l)'^f. Notice the similarity with 
the also ring-like distribution of vevs. In particular, an infinitely thin ring aX p = pq will 
be considered, although the construction can be easily generahzed to other configurations 
with the same symmetry. The choice of masses fixes the distribution of branes along the 
p, 02 directions. We still have to fix their distribution in 6', (p. Since in the low energy theory 
these coordinates play no role, one can just consider the simplest possibility, i.e., a uniform 
distribution such that the S'0(3) symmetry is also kept. This dramatically simplifies the 
task of writing an ansatz. Figure [T] depicts the situation. 




Figure 1: A scheme of the configuration on the p, 02 plane, which is the complex -u-plane up 
to a factor of A. The inner circle at radius po is the distribution of Nc vevs for the adjoint 
scalars. The outer circle represents the Nf masses of the fundamental hypermultiplets. 
The associated flavor branes are also extended along the transverse a, 0i directions. These 
smeared distributions restore the U{1)r broken by each single vev or mass. The point at a 
scale Pq represents a flavor probe brane whose dynamics will be studied in section 0J 

The D5 flavor branes are magnetic sources for the F(3) RR-fleld strength. As in |12j . 
a supersymmetric solution can be found by modifying the Bianchi identity for the F(3) 
and inserting the suitable ansatz in the supersymmetry variation equations ()2.12|) . The 
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worldvolume action for the flavor branes is: 



s 



flavor 



d xe^ 



m + / ^[^6)] . 



(3.2) 



We are interested in the WZ part, which, being an infinite sum can be promoted to a ten- 
dimensional integral. 



(3.3) 



where Vol{y^ is the brane density, subject to the normalization condition J Vol{y^ = Nf. 
For the distribution described above: 



^0/(3^4) 



8^ 



^(P ~ Pq) sin 9dp A d(f)2 A d9 Ad(p 



The modified Bianchi identity reads^: 



Nf 



dF(^s) = 2Kl^)%Vol{y^) = g,a' -^6{p 



sin 9dp A d(j)2 Ad9 A dip 



(3.4) 



(3.5) 



where ()2.7|) and ()3.4p were used to obtain the second equality. The natural ansatz satisfying 
this condition that generalizes (j2.1(jp is: 



'(3) 



Ncgsa' 



-9 



A dp A {d(f)i + cos 9dip) — gd(f)2 Ada A {dcpi + cos 9dip) + 



Nf ^ ^ ■ 

+ {g + - Pq)) sm9d^2 Ad9 Ad(p 



(3.6) 



where G is the Heaviside step function. Maintaining the same ansatz for the metric^ ()2.9p 
and the expressions for the Killing spinors ()2.14|l . ()2.15j) . one finds a very slight modification 

"'in general, if for a form = c?A(„_i') there is an action / \/\g'W^ + / G A A, the equation of 
motion for the form reads d * F — sign{g){—l)^^"^^G. In this case, the relevant part of the action (go 
to string frame for this computation) reads: — ^Ji^ — ^ J -y/jglF^^^ + T5 J 1^0/(3^4) A C(g) so the equation of 

motion is ^^l^ d * F(j^ = —T5 1/0^(3^4). Taking into account F(3) — — * Fjy) and we arrive at H3.5(l . 

'If one considers a more general form for the metric 

1 



-dxf r^ + z{de^ + sin^ 6*^^^)+ 



gsNcW 



We-^'^(dp^ + p^dcfyl) + (3^ ( da^ + a'^(d(j)i + cosddif) 

z \ 



where a{p^ a), (3{pt a), the Killing spinor equations impose that both a and /3 are constant. 
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of (ITTT)|) - (jT^ . namely: 

El 

2N, 



Nf 

9 + = -P^\ (3.7) 



e^* = -, (3i 



zz 

g' = -2e-2*p(T<l) , (3.9) 
'"^'"Mg + TTW^iP - Pq)) + 2;.-Vae-2*$' . (3.10) 



-z 



2N, 



c 



As in the unfiavored case, ()3.10p is implied by ()3.7|) . ()3.8p and it is possible to check that these 
equations ensure the equation of motion for the 3-form (i(e* *-F(3)) = 0. The generalization 
of T27K^ is 

Nf 

""WTF^iP - Pq) + - = ^(P^' + ^' + P^") • (3-11) 



' c 



This equation seems very difficult, if possible at all, to solve in general. However, the 
important point is that it uniquely defines the relevant solution. For p < pg, the solution 
should be equal to the unfiavored one, which is explicitly known 1)2.211) . A first argument to 
see this comes from an electrostatic analogy: a spherically symmetric distribution of charges 
does not alter the field in its interior. Such a reasoning was used in j2ZI in order to find 
the supergravity solution corresponding to shells of five-branes uniformly distributed on 
spheres. Notice the similarity between such a setup and the one discussed in this note. A 
second argument comes from field theory and holomorphic decoupling (see the final part of 
section 1221) • 

Thus, we are left with the problem of finding the function z{p,a) for p > pq. It will be 
shown in section that z(p, 0) can be written in a simple form and leads to several gauge 
theory predictions. By continuity, and regarding the previous paragraph z{pQ,a) must be 
the same as in the unfiavored case and therefore is known. From ()3.7)1 . one can read how the 
derivative changes when traversing the brane shell, fixing z'{pQ,a). This set of boundary 
conditions fixes the solution. In figure 121 a numerical approximation to the function z{p,a) 
in two particular cases is presented. 

There is a last subtle point that requires discussion before closing this section. We have 
introduced a number of fundamentals and found new equations that define a new solution. If 
this new solution still describes a theory with A*",, colors, equation ()2.5|) should still hold. Let 
us now prove so. Comparing to ()2.(ij) . the integrals in 0i, 02 are trivial, but the angle 6 is not 
manifest in the geometry any more. This problem is solved by noticing that the components 
of -F(3) relevant to this integration (those not containing d6 or dip can be written as (see 
()3.(i|l ) NcQsCt'dg A d(j)2 A dcpi and that, at infinity, the limits of integration in 6 correspond to 
(cr = 0, p = oo) and {a = oo, p = 0). Then: 

J F(3) = 47r^5(,A'ca'(5(|((<^=oo,p=o) - g\{a=o,p=oo)) = ^Tv^gsNca' (3.12) 
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Figure 2: Numerical approximations to the function z{p^ a) that enters the metric. Both 
plotted solutions correspond to fixing c = (po = V^) and pg = 5. It is apparent from the 
graphics the qualitative change of behavior at such values of p. On the left Nf = Nc. On 
the right, Nf = QN^, a non-asymptotically free field theory that eventually hits a Landau 
pole. This corresponds to z becoming zero at some value of p, as shown in the graph. 



where the equalities (y(|(o.=oo,p=o) = and (7|(o-=o,p=oo) = — 1 have been used, and also the fact 
that g is continuous along the integration path. The first equality is trivial since for p < pq 
the solution coincides with the unfiavored one so g can be read from ()2.22j) setting ^ = 0. 
The second one will be proved in section ESI Thus, since ()3.12p coincides with ()2.6|) . the 
quantization condition still holds in the fiavored case. 



3.3 Gauge theory features 

Let us now study the gauge theory implications of this computation. The p < po region is 
unchanged with respect to section |21 when the fundamentals are present, so in the following, 
only the p > Po region is considered. By repeating the argument of [201 reviewed in section 
123 it is straightforward to generalize (j2I27j), (IT^ to: 



9ym 



9 



Ma=o), (3.13) 
Nf 



YM 



2(^?U=o + ^e(p - pg))iV,02 (3.14) 



Equation ()3.9|) implies that g\a=o is a constant. Repeating the argument that the solution at 
small p should be the same as the unfiavored one, we find (as in ()2.24|) ) g\a=o = — 1- Thus: 

9yM = -2iV,(l - ^e(p-pQ))02 (3.15) 
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One can compute the monodromy of the 9ym angle when making a full U{1)r rotation 
through the space of vevs (at a fixed modulus pA): 

02 -.02 + 477 ^ ^^^_2(2iV,-iV;e(p-pQ)) (3.16) 

which depends on NfQ{p — pq), the number of fundamental hypermultiplet masses encircled 
by the path. This result is in exact agreement with the field theory expectations (for a 
review, see |2H1)- 

Let us now turn to gvM- Using g\a=o = —1, it is immediate to integrate ()3.13|1 . Imposing 
that z matches the unflavored solution for p < pq and is continuous at p = pg, one obtains: 



(iv, - ^e(p - pq)) logp + ^e(p - pq) logpe] (3.17) 



91m 47r2 

Using ()2.29|1 . the /3-function at an energy scale p is straightforwardly computed: P{gYM) (p) = 
— %^{Nc — ) where Nf{p) is defined as the number of flavors for which the modulus 
of their masses is smaller than the scale. The fact that matter fields with bigger mass do 
not contribute is usually an approximate statement, but, due to the radial symmetry, in 
this case it is exact as will be shown below. Their effect is just to modify the dynamically 
generated scale A. For Nf > 2Nc, the theory loses asymptotic freedom and for Nf > 2Nc, 
it eventually hits a Landau pole. In the gravity solution, this is translated in the fact that 
when Nf > 2Nc, the function z vanishes at some finite p and the metric becomes singular. 
An analogous behavior in a D3D7 solution was discussed in jllj . 

The holomorphic coupling can be written by using ()3.15j) and ()3.17|) (define Uq = PqA): 

r = t — log — , for m < Mo , 

TC A 

T = I — log — , for Mo < |m| < Mq , 

TC A 

r = z-^log-^ + -(Ar,--^)log-, for|n|>nQ, (3.18) 
/vr A vr 2 A 

This behavior precisely matches the field theory result. The proof is a simple generalization 
of the one of ^H] for the unflavored case. The coupling is given in terms of the so-called 
prepotential. As argued in [21], non-perturbative contributions to the prepotential vanish in 
the large N limit as long as one probes the theory far enough from the vevs |m — aj| ^ N~^ 
(or, in this case, also from the masses of the hypermultiplets). The perturbative prepotential 
can be read, for instance, from Slightly adapting conventions in order to match those 

of HE!: 

^ = ^ - o.)^ - ^ 5: - "")^ ^^4?^ (3.19) 
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The gauge group is broken to U^l)^" due to the vevs of the scalars. We now probe the 
theory by considering all but one of the vevs in the ring-like distribution and moving around 
the other one, u, so we have [7(1)^'="^ x f/(l) where the last U{1) corresponds to the probe 
brane. In the large Nc limit, the coupling is: 

92 J" i ^ (n-ai)2 i ^ (^u-rriay 
rl-u = — — = — > log^ — -TT^ > log^ — — — (3.20 

i=l a=l 

In the large N^, Nf limit, the sum over discrete distributions can be replaced by integrals, 
namely: 

r{u) = ^ / " ap<*>(a)log — -j^ 4^ / ^PrriQim) log — (3.21) 

where p<*>(a), pmQijn) are the density distributions of vevs for the adjoints and masses for 
the fundamentals respectively. The densities must satisfy the normalization j (Pap^x^,^ (a) = 
Nc-, j (PmpmQ{m) = Nf. Therefore, the ring-like distributions are p<*>(a) = -^^Sda] —uq) 
and pmoim) = 2^^5(|m| — uq). Performing the integrals ()3.21|1 and adjusting an additive 

Q ^^^^^^ 

constant that can be reabsorbed as a rescaling of A one finds precisely ()3.18|) . The rescahng 
of A is what one usually obtains from holomorphic decoupling. 

In the context of fractional branes at orbifolds with massless fundamental hypermultiplets, 
results similar to those reported in this section were presented in |S] . 

4 Mesonic excitations 

This section is devoted to the analysis of the meson-like excitations, both in the quenched 
and unquenched backgrounds. In order to address this question, let us add an additional 
flavor to the theory with a mass defined by its position in p, 02- The flavor group is now 
t/(l)'^/+i. Obviously, the effect in the geometry of this new brane is suppressed by a power 
A^~^ and is negligible. Thus, its excitations can be described by considering it a probe in 
the background described in the previous sections. 

The spectrum can be computed by expanding the action to quadratic order in the fluc- 
tuations. The analysis is similar to that in (Hj, where also an A/" = 2 theory was analyzed. 
There are, however, important differences in the low energy field theory: apart from having 
a large number of hypermultiplets, here there is no massless adjoint hypermultiplet. 

Let us consider a probe brane extended along xq, Xi, X2, Xs, cr, 0i and sitting at a point 
defined by 9g,ipg, pq,(f)2q- Its worldvolume action reads (below, string frame will be used, so 
the metric ()2.9j) should be multiplied by et): 

Sd5 = tJ- f d'^e-''^-P[g] + 2na'F+ j P[C(6)] + j P[C(2)] ^ F ^ f\ (4.1) 
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although the last term does not contribute at quadratic order, as can be seen by choosing 



a gauge in which C(2) = QsNcd' [g + ^6(p — pq) j (i02 A (d^i + cosOdcp). We need the 



expression for the C(6) in the background which is defined as (iC(6) 
the Hodge dual of ()3.6|) . one finds: 



*F(3). Taking 



F(7) = gsNca'dxo A (ixi A dx2 A (ixa A 



/ 4<1? 

^ sin eirffii Ad(fAda + 

pa 



gz-e 



in 9d9 A dip A dp - [g + ^^Q{p - Pq)) -^dp Ada A {d(j)i + cos 9d(p) (4.2 



pa 



sm 



2N, 



pz^ 



Using ()3.7|) - ()3.10|) . one can write the following simple expression for the RR-potential: 



C{Q) = gsNca'dxQ A dxi A dx2 A dx^ A ( ze sin 6d6 A dip da A {d(j)i + cos 9dip) ) (4.3) 



4.1 Fluctuation of the scalar fields 

In order to describe the small fiuctuations, the simplest method is to consider the static 
gauge and use x^, a, 0i as worldvolume coordinates. The embedding can be written as: 



9 = 9g + 66 , ip = iPg + 5ip , P = pg+ 5p 



+ 



(4.4) 



where all the fiuctuations depend on x^, a, 0i. It is immediate to see that at quadratic order, 
these scalar fiuctuations do not mix with those of the worldvolume gauge field which will be 
studied in the next section. Expanding the action 1)4.11) up to quadratic order one finds: 

/r X X X 

d^xd<l),da --e-''''g,N,a\dJpf + -a{d^5pf + —{d^,5pf + 
12, z I la 



1 oO" 



-^-A-e-^^g.^c(J{dM2Y + -JAdMif + — p: 



la' 



1 ~ 1 ~ 7P^ 

^-g^N,ola{dJQf + -az^^{djdf + ^{d^M'' 
I I la 

1 _ 1 _ rp- 

^-g,N,ola^xY?Qg{dJp^f + -aze^^ ^mHg{dJp^f + sin2e,(9^,5<^)2 + 
I I la 



,2$ 



(4.5) 



where z, $ should be understood as functions of cr at p = pg. We have defined [d^SpY = 
r]^^{d^5p){du5p) where p^v span the four Minkowski directions. The fiuctuations of 6p,6(j)2 
are already decoupled and it is very easy to decouple those for 66,6ip. Let us expand the 
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fields as: 



^02 = X2(f^)e*^'''cos(/0i) , 

66 = ^±^^1±^ sin 9,e^'-Ml4>,), 

5^ = ^±^^l^^e^'-oo.m. (4.6) 



The equation of motion for xi, X2 reads: 



i2 

M\, + ^,{a^^X^)--X^ = ^ ■ (4.7) 



ze2* a 

wliere we liave defined tfie four-dimensional mass = —k^k^ and = QgNca'M'^. 
The equations for x± read: 

M\± + -dAaze"'d^x±) - -^X± T -dAze''')x± = (4i 

4.2 Fluctuation of the gauge field 

Tlie equation of motion for the gauge field at linear order is: 



da{e-''^-det{P[g])F-') = (4.9) 

where a, b run over the six worldvolume coordinates and, again, /x, u over the Minkowski 
four-space and are contracted with the flat metric 77^,^. Let us now choose a Lorentz gauge 
df^A'^ = 0. The following relation is found: 

d^,A^, = -ad^iaA^) (4.10) 

Defining a transverse polarization tensor and expanding: 

= e^X3(f^)e^^''^cos(/0i) , 

= -X4(a)e^'=-"cos(/0i) . (4.11) 
a 

one finds from ()4.9|) that Xa^Xi also satisfy equation ()4.7p . 
4.3 Analysis of the spectrum 

It is interesting to notice certain similarity of the equations ()4.7p . ()4.8p with (3.6) and (3.31) 

of P , respectively. Thus, the analysis of the spectrum is also similar to the one in P . 
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However, the / we have used here is not the eigenvalue of some operator acting on the 
SU{2)r generators, but it is the charge under a U{l)j C SU{2)r, i.e., the eigenvalue of a J3 
operator. The relation between the two would require a better understanding. The bosonic 
content of the generic {I > 2) Af = 2 massive vector multiplet is given by a massive vector 
and three real scalars in the | of SU{2)ji and two real scalars in the ^ and 

Since fj4.7|) is satisfied by Xi, "^ = 1,2, 3, 4, we have in fact a massive vector and three real 
scalars with the same mass in the same representation of SU{2)ji. In order to complete the 
multiplet, two more real scalars with the same mass are needed. They are the x± of ()4.8|) . 
This can be shown by mapping equation ()4.7|) to ()4.8|) . following a procedure used in jSJ for 
a similar case. In fact: 

X+,{i=L) = (r'^'^d^{xi,ii=L+i)(^^^^) 
X-,ii=L) = t^^"'5.(x.,a=L-i)^'"^) (4.12) 

In order to prove this, substitute in (j4.7|) I ^ L ±1 and define F± = a^^^Xi- Multiply the 
equation by ze^^a^^ and derive with respect to a. Now, substituting d^F^ = x±^^^^ and 
multiplying by a^^^^ , one arrives at ()4.8|1 with / = L. 

Thus, the bosonic fluctuations yield the bosonic content of a tower of A/" = 2 massive 
vector multiplets which can be identified with the mesonic excitations of the theory. The 
fermionic spectrum can be inferred from supersymmetry or could be directly computed along 
the lines of [32] • 

The physical, mesonic, excitations should be those satisfying equations ()4.7j) or ()4.8j) and 
moreover being regular at cr = and normalizable as a — 00. Normalizability is a problem- 
atic issue in this setup due to ill-behaved UV and will be discussed in section 1^711 

About regularity at the origin, one would need that the x functions are finite and even 
(odd) around a = when I is even (odd). From ()4.7|1 . one sees that near a = 0, Xi,i = 
cicx' -|- C2cr^'. Clearly, regularity selects C2 = and Xi,i which in fact satisfies the 

criterion of parity stated above. For the x+ modes, substituting in ()4.12|) . we find x+,i ^ 
which is also fine with the parity condition. Checking the X- modes, involves an extra 
subtlety since substituting the leading Xi behavior in ()4.12j) one just gets zero. But from 
()4.7|) . noting that near a = 0, i oc a we see that both z and e^* behave as ki + ^20"^ where 
ki, ^2 are some constants. Thus, the subleading behavior of Xi is described by ci(j'(l + fcscr^), 
so substituting in ()4.12j) one also finds X-,i oc cr' which, again, is fine. 

4.4 Estimate of the mass spectrum 

In D5-brane solutions, like the one we are dealing with or the Maldacena-Nunez background, 
the dilaton diverges in the UV. The supergravity formalism cannot be trusted in that region, 
and the UV completion of the field theory is a little string theory. Therefore, it is not a 
surprise that there are problems in the UV when trying to determine the spectrum. In 
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fact, there is no normalizable mode. However, physically, one would expect to have some 
modification of the UV that cures this problem, since in the dual theory there are, of course, 
physical excitations. In this section, a regularization procedure similar to those in Oil 
proposed, which basically consists of cutting off the ill-behaved region. It is curious to notice 
the qualitative difference with where glueballs in Maldacena-Niiiiez were considered. 
In that case, even though there was no normalizable mode, the authors found that there 
existed one leading and one subleading mode in the UV and identified the subleading as the 
physical one (in a different context, the same criterion was also proposed in [SE])- In the 
present case, there is no leading UV mode because the two modes behave as sine and cosine 
of some function of a. 

In order to justify the procedure, it is interesting to notice the similarity in the UV of this 
setup to the D5D5 brane intersection in flat space. Of course, in such a case, there are also 
no normalizable modes p?T| . But in the UV, by a chain of dualities one can connect the 
system to an M5M5 intersection which indeed has physical normalizable oscillations j|37j. 
Let us deflne: 

y = \oga (4.13) 
such that equation ()4.7p can be written in the Schrodinger form: 



dy 

The potential reads: 



% - V{y)x^ = . (4.14) 



e 



2?/ 



V{y) =1'- M'—^ (4.15) 



where y G (— oo, oo). The function ze^"^ is implicitly known in the quenched case (see section 
12)) and can be determined numerically in the unquenched described in 13.21 

For D5D5 intersecting in flat space [37j the potential reads: 



V{y) = f-M'-^^ (4.16) 

Figure El shows the qualitative similarity between the two cases. For simplicity, the potential 
()4.15|1 is plotted for the quenched case, although the qualitative behavior for the unquenched 
one is the same. From the flgure, we see that starting from y k, A ^ a ^ 50, all the potentials 
coincide, i.e., the IR effects are no longer important. It is natural to choose such a scale as a 
cutoff. In the plot, the UV problem is apparent since Viy oo) < and such a Schrodinger 
potential cannot have a discrete spectrum. 

The masses can be estimated by applying the WKB method, which should be reliable for 
large excitation number. Let us focus on equation ()4.7|) . In the notation of appendix El we 
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m 



the left represents the potential (j4.15|) with fixed 
1 for different values of pq = 1,2,3,4 (the biggest the pg 



Figure 3: The figure 
^ = 0, 1 = 2, M = 3, c 
the upper is the line). The thick line corresponds to ()4.16p . The figure on the right is for 
fixed ^ = 0, 1 = 2, M = 3, Pg = 3forc = — 1,0,1. The different lines are almost coincident 
and the meson masses should not strongly depend on c. 



have: 



a 



ze 



2'I> 



P=-- 



a 



h = h 

S2 = 1 , 
V\ = P2 



= Si 

r2 -- 



-- ri 
-1 



-I' 



S3 = rs 



(4.17) 



where it was used that 2;e^* ~ (up to logarithms) for large a. In fact, the problem for 
large a appears in this formalism due to the fact that ri — r2 — 2 is zero, but it should be 
a positive number in order to apply the WKB. Now, suppose that we can change the UV 
by that of the M5M5 intersection, which would have r2 = —2 37^ and leave all the rest 
unchanged. With this, the WKB approximation to the masses reads ()A.5|) : 



TT 



WKB ~ -^n{n - 1 + 3/) , n > 1 



where flA.6 




cutoff 



(4.18) 



(4.19) 



Now, the problem for large a appears in the fact that this integral is divergent if (Tcutoff ^ oo. 
As said above, a natural scale to cut the integral is around y ^ 4 ^ a ^ 50. Although this 
involves some arbitrariness, everything is fixed in terms of a single parameter (J cutoff and thus, 
it should be possible to estimate how the meson masses depend on the physical parameters 

C, pq,Nf,PQ. 
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Figure m shows such a dependence in the quenched case, which is easier to deal with since 
there is an expression for z ()2.2H1 and no numerical integration of a system of PDEs is 
needed. The quantity proportional to the meson mass, is plotted versus pq. The graph 
on the left shows how changing (Jcutoff does not change the qualitative properties of the 
graph, although of course shifts ^. The graph on the right, shows the behavior for different 
values of c (related to po by ()2.23p ). For ^ Po; the curves become degenerate as expected 
on general grounds. The minimal value for the meson masses is attained at Pq = po. It is 
natural to conjecture that this is due to the fact that when < ^ >= ruq for some of the 
eigenvalues, there are cancellations between the superpotential term ()3.1|) and the mass term 
yielding effectively massless quarks. 




Figure 4: On the left, ^ ^ versus Pg for c = 0. Starting from above, the three lines correspond 
to (Tcutoff = 50, 70, 90. On the right, versus Pq with fixed Ccutoff = 50 for c = 0, 2, 4. In 
both plots, — > 0. 

An important question is how the quantum effects produced by unquenched flavor affect 
the spectrum of physical excitations. From the discussion of section 13.21 it is clear that for 
Pq < Pq the spectrum is unchanged by the unquenched fundamentals. When pg > pq, the 
computation requires a delicate numerical analysis of the PDF p.llj) in order to obtain the 
function 2;e^*((j) for given pq, pg, c). I could not find any significative change of the 
meson masses when the Nf ~ hypers are introduced, even when pq > pq. However, this 
conclusion should be taken with a grain of salt due to the difficulty of obtaining a precise 
numerical solution near a = and also for pg^ pq. 
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A Appendix: WKB approximation, useful formulae 

This appendix collects the results and notation of concerning the WKB approximation. 
They have been used in section 14. 4[ Consider a differential equation: 

dMi^)da<P) + {M'h{a) + p{a))<p = (A.l) 

where the functions f,h,p behave near a ^ 0, a ^ cxd as: 

/ ^ fia'' , h ^ hia'^ , p ^ picr'' , (fT ^ 0) 
f^f2a'\ h^h2a'\ p^P2<y'\ (a^oo). (A.2) 

The WKB approximation is consistent provided S2 — si -|- 2; ri — r2 — 2 are positive numbers 
and S3 — Si -|- 2; ri — r3 — 2 are positive or zero. Define: 

ai = S2-Si + 2, /3i=ri-r2-2 (A.3) 

and 



a2 = |si-l| or = ^/(si - 1)^ - 4— (if S3 - Si 2 = 0) 



/32 = |ri-l| or /32 = J(ri-l)2-4^ (if n - r3 - 2 = 0) (A.4) 



The masses are approximated by: 



= ^n{n - 1 + ^ + ^) + 0(^0) , n>l (A.5) 



with: 



00 







da J J (A.6) 
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